Brownian motor in a granular medium 



R. Balzan/'2 F. Dalton,^ V. Loreto/'^ A. Petri,^ and G. Pontuale^ 

^ Dipartimento di Fisica, Sapienza Universita di Roma, 
Piazzale A. Moro 2, 00185 Roma, Italy 
^Consiglio Nazionale delle Ricerche, 
Istituto dei Sistemi Complessi, sede di Tor Vergata, 
Via del Fosso del Cavaliere 100, 00133 Roma, Italy 
^ISI Foundation, Viale Settimio Severn 65, 
Villa Gualino, 1-10133 Torino, Italy 

Abstract 

In this work we experimentally study the behavior of a freely-rotating asymmetric probe immersed 
in a vibrated granular medium. For a wide variety of vibration conditions the probe exhibits 
a steady rotation whose direction is constant with respect to the asymmetry. By changing the 
vibration amplitude and by filtering the noise in different frequency bands we show that the velocity 
of rotation does not depend only on the RMS acceleration F, but also on the amount of energy 
provided to two separate frequency bands which are revealed to be important for the dynamics of 
the granular medium: the first band governs the transfer of energy from the grains to the probe, 
and the second affects the dynamics by altering the viscosity of the vibro-fluidized material. 

PACS numbers: 45.70.-n, 05.40.-a, 05.70.Ln 
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Exploiting spontaneous motion for producing energy has been a long standing dream, the 
impossibility of which in equilibrium systems was definitively established with the advent 
of thermodynamics showing that perpetuum mobile of the second kind is contrary to the 
second principle for any closed system. At the microscopic level this conclusion relies on 
molecular chaos as discussed by Smoluchowski and Feynman , though the laws of 
thermodynamics do not prevent the possibility of extracting work from systems that are 
isothermic and in a stationary state but not at equilibrium. Such a state would require some 
sort of spontaneous symmetry breaking and constitutes the subject of Brownian motors 
Experimental realizations, generally named thermal ratchets Jiave been recently studied in 
different situations, ranging from the nano to the micro scale 3|jj4|. Macroscopic realizations 
have been developed in the field of granular media (GM) js, 6|, a field of great relevance 
both in applied science, for their industrial relevance, and in theory, for their challenging 
properties and behavior ^\. Strictly speaking granular systems are athermal since the motion 
of the elementary constituents of the medium, the grains, is not affected by the ambient 
temperature. It may therefore seem contradictory to speak of thermal ratchets; however 
grains may gain kinetic energy by external mechanical perturbations, and so be ascribed a 
'temperature' defined by their motion j^, yet never reach equilibrium because of the presence 
of inelastic collisions and friction. Therefore, GMs can be used to implement macroscopic 
.ealizationscf thecal .a.che.s. both expe..e„ta„y flBQ and theoretically UU^ 

The experimental realizations cited above were principally targeted towards the observa- 
tion of spontaneous collective oriented motion of the grains, whereas in this work we wish 
to focus on the spontaneous motion of an external asymmetric intruder in the absence of 
any collective granular motion. Numerical simulations have demonstrated that such a phe- 



nomenon is possible 



We present novel results showing that the chaotic granular motion 



can indeed propel the asymmetric probe in a persistent direction. More specifically, the 
immersed asymmetric probe, under suitable conditions of fiuidization and viscosity, exhibits 
consistent and significant motion in a constant direction with respect to its asymmetry. 
Under different conditions the system can exhibit collective motion of the grains (e.g. con- 
vection) ; however in these cases the symmetry breaking is performed by the GM and so the 
direction of asymmetry of the probe is irrelevant. 

The main components of the experimental apparatus are the granular medium, the shaker 
and the probe as shown in Fig.[Il The signal to the shaker is supplied by a function generator 
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Figure 1: Example of the symmetric and asymmetric probes used (a). The probe is attached to a 
freely-rotating support (b) and is inserted into a beaker (c) containing the beads. The beaker is 
vibrated by a shaker (d) (Bruel & Kjaer type 4809) and the instantaneous angular position of the 
probe is registered by a rotary encoder (e) with resolution 1/500°. 

and is amplified and filtered with an efficiency of 22 — 23 db/oct. The vibration was applied 
in the range from 10 to 1000 Hz and the system acceleration (0 < F < 4, normalized to 
gravity g) is measured by an accelerometer fixed to the vibrating stage of the shaker. The 
probes, each 3 cm tall, are immersed to a depth H in the GM (glass beads of diameter 
2 mm ± 10 %) which approximately half-fills a beaker, 90 mm diameter and 120 mm tall. 
It is well known that when a GM is shaken with a periodic signal, a wide range of coherent 

• n 

motions can arise [16||. Generally speaking, such motions are due to a dynamic symmetry 
breaking in the system Q| and vary from the presence of isolated and discontinuous collective 
movements of a minority of beads to major effects like inclination of the surface or formation 
of a single vortex in which the entire GM rotates. In our system these coherent motions 
typically arise when the shaker is excited by a pure square or sine wave in a range of 
frequencies up to about 150 — 200 if 2;. For example, well defined internal convective motion 
was observed which caused the probe to rotate until it stops at a specific angle, irrespective 
of the starting point. Alternatively, in the presence of vortices the probe was seen to rotate 
indefinitely in a constant direction at roughly a constant velocity, with no dependence on 
either the presence or the direction of the probe asymmetry - the probe clearly follows the 
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macroscopic collective motion of the granular material. However, above 200 Hz or with 
white noise excitation, only chaotic motion of beads was observed with no collective effect, 
while over 1000 Hz the granular material seems no longer to respond to the agitation even 
for r ^ 1. Hence white noise was applied to the shaker in order to eliminate any convective 
or collective motion of the GM, and this allows us to apply energy also below 200 Hz without 
triggering collective motion. The use of non-periodic signals requires a correct definition of 



in which z (t) is the instantaneous acceleration of the granular container and (...) is a time 
average. T is also expressed as a function of the Power Spectral Density (PSD) where ( (u) is 
the Fourier transform of z{t). Equation ([T]) reduces to the usual definition of F for sinusoidal 
excitations. Other definitions also exist which reduce to the usual F for periodic signals, 
however the above choice is motivated by the fact that Eq. ^ is proportional to the overall 
energy provided to the granular medium, and furthermore that it yields the value F = 1 at 
an evident transition point in our experiments (this transition is also observed at a calculated 
value of F = 1 with pure sinusoidal excitation). 

Using white noise in a large enough frequency range eliminates convective motions and 
the granular exhibits a uniform and homogeneous behavior. Under these conditions it may 
be expected that no net rotation of the probe would be observed due to the random nature 
of the action of the granular material on the probe surface. Contrarily, however, there is 
an almost regular rotation of the probe when F > 1, the direction of which depends on 
the orientation of the probe asymmetry. No net rotation was observed using the symmetric 
probes. The response of the probe at various excitation intensities and frequency ranges has 
been studied with the observation that the probe velocity is dependent on the upper and 
lower frequency limits of the shaking signal. The probe immersion depth H into the GM 
has also been varied. 

Figure [2] presents a general result illustrating the asymmetry-dependent mean rotation 
velocity {u) of the probe as a function of the vibration intensity F. Below F = 1 the velocity 
{uj) = 0. There is a clear transition at F = 1 at which point uj initially increases rapidly 
before reaching an asymptotic value which is dependent on the probe used (e.g. probes with 
4 teeth rotate faster than those with 2), and the direction of its asymmetry - inverting 
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Figure 2: Probe average velocity as a function of the shaking acceleration for the 2 and 4 teeth 
asymmetric probes (numbered 1 to 4 as in fig. [I]). The vibration signal is white noise filtered in the 
range 15 — 300 Hz. Inset: Average rotation velocity u; as a function of the immersion depth H of 
the probe at T = 2. At H = 3 cm the probe is fully immersed in the medium. 

the probe asymmetry results in a simple inversion of u. Notably, symmetric probes show 

(uj) = ovr. 

In the inset to Fig. [2] the average velocity of the probe for various immersion depths H 
of the probe into the granular material is shown. For small H < 1 cm there is no rotation 
and we assume that the action of the GM on the probe is insufficient to overcome friction 
in the support. 

Increasing H causes the probe surface to engage more completely with the GM, and the 
velocity increases roughly linearly until H ^ 2.5 cm. At this point the probe is almost 
fully immersed in the medium and the velocity reaches a brief plateau. Further immersion 
H > 3 cm causes the rotation velocity to fall quickly to 0. This behavior suggests that the 
region of the GM which applies the spontaneous torque to the probe is close to the surface. 

By changing the intensity and the frequency range of the applied white noise it may be 
seen that the behavior of the system is principally determined by the energy provided in two 
major frequency bands. We consider a first series of experiments in which we apply a single 
frequency band to the GM by applying a band-pass filter to white noise from /i to /2 Hz. 
Firstly by altering 15 < /i < 60 Hz with fixed /2 = 300 Hz then by fixing /i = 15 Hz and 
altering 300 < /2 < 900 Hz. As the frequency band is widened, the overall energy provided 
to the system would increase; to investigate therefore only the change in frequency, the signal 
intensity is reduced in order to maintain a constant overall energy input, proportional to T 



(the probe behaves similar to Fig. [2] as a function of F). 

For different /i and fixed F the probe approaches a different asymptotic value. Figure [3] 
displays this asymptotic value of the mean probe angular velocity {u) (taken at F = 4) 
as a function of /i, the lower band-pass cutoff. It is clear that the {u) decreases slowly 
from /i = 15 Hz to /i = 30 Hz, from 30 to 40 Hz there is a jump which corresponds to 
an absorption peak in the apparatus response function [2^, due to the granular material 
efficiently absorbing energy from the vibration and so mobilizing the probe. For /i > 
40 Hz (uj) continues to decrease slowly. Importantly, this jump between 30 and 40 Hz is 
independent of the value of the higher cut-off frequency /2 (in this example fixed at 300 Hz, 
though fixed at up to 900 Hz is other experiments). Thus we conclude that the frequency 
band from roughly 30 to 40 Hz is crucial to generating motion of the probe. 

We now fix /i = 15 Hz and increase the upper band-pass cutoff frequency 300 < f-2 < 
^00 Hz, still maintaining a constant overall energy input by reducing the signal amplitude as 
necessary, as schematized in the inset to Fig. [4](a). The main Fig. [4] demonstrates that the 
mean velocity actually decreases as the frequency band is widened (black circles). This may 
seem contrary to Fig. [2] but in fact is a direct consequence of Eq. ^ and the reduction in 
signal intensity, specifically in the 30 — AO Hz band, which is required to maintain a constant 
F. In fact, if one divides the velocity obtained by the amplitude applied (relative to that 
with /2 = 300 Hz), one obtains an almost constant velocity (green triangles), indicating 
linearity of uj with 7. 

To overcome this effect, we consider a final series of experiments in which the amplitude 
of the signal sent to the shaker is kept constant. In order to supply sufficient energy in 
the bands which interested us (15 — 60 Hz and 200 — 900 Hz), it was necessary to apply 
a band-reject from 60 to 200 Hz; the spectrum applied to the shaker is schematized in the 
inset to Fig. [4](b). In this manner a constant energy input in the lower band is maintained 
while the upper band can be independently widened. 

In Fig. [4] we show the variation of the asymptotic probe velocity as the upper band is 
widened from 300 to QOO Hz (blue x). The velocity is seen to increase only slightly despite 
the large increase in the band width. 

When taken together with the observation that motion ceases when energy is removed 
from the lower-band irrespective of the upper band, then the logical conclusion is that the 
upper frequency band is not responsible for the spontaneous torque causing the rotation. 



6 



A 0.15 



Spectrum (schematic) 



— t5Hz 

— 30Hz 

— 60Hz 



100 200 

Frequency (Hz) 



30 40 50 

Lower cutoff frequency (Hz) 




Figure 3: The asymptotic value of the probe velocity (i.e., (uj) at F = 4) while varying the lower 
cutoff /i of the band-pass filter from 15 to 60 Hz, necessitating an increase in the filter amplitude 
in order to maintain a constant F. The mean velocity decreases approximately proportional to fi 
with a step between 30 and 40 Hz, coincident with the absorption peak in the response function. 

The slight increase in velocity, therefore, may be ascribed to a reduction in the viscosity 
of the GM; indeed we imagine a high-frequency vibration merely vv^eakening the contacts 
between grains which, however, remain caged in place by their neighbors, whereas a low- 
frequency motion actually imparts macroscopic motion to the grains. Such a phenomenon 
has already been observed in [8|. 

The viscous torque Mg can be written as the product of the probe viscous coefficient in 
the GM ub and the probe velocity u: 



-UbUJ. 



(2) 



where ub is explicitly dependent on the excitation band, and in fact is inversely proportional 
to r _ 

aiPSDif)) 



r 



where a{PSD{f)) is some function of the spectrum applied to the shaker. 

Thus we infer that, when the lower cut-off is increased, F is drastically reduced, which 
both increases the viscosity, and reduces the spontaneous torque F generated by the GM. 
When the upper cut-off is increased, however, F remains almost constant but a decrease in 
viscosity is observed due to the presence of the term a{PSD{f)) in Eq. (JSj). 

In summary, the origin of the observed behavior can be traced back to both the mechanical 
absorption of the GM and its fluidization. In the case of F < 1 rotation does not occur as the 
GM is insufficiently mobilized and maintains a high viscous coefficient ub- The reduction 
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Figure 4: The value of the probe velocity (measured at F = 4) while varying the upper cutoff /2 of 
the band pass filter. This was performed both with a single band, maintaining the overall signal 
energy by decreasing the amplitude as appropriate, and with a double-band at constant amplitude 
(see inset). In the first case, the velocity decreases (black circles) displaying a linear behavior with 
respect to the amplitude (green triangles). In the second case, increasing 7- fold the energy in the 
upper band produces only a slight increase in the mean probe velocity (blue x); this is in contrast 
to the result in Fig. [3] where a reduction of approx. 15% resulted in a 75% reduction in the velocity. 

of ub with r > 1 confers sufficient mobility to the probe and so rotation occurs. It is clear 
from Fig. [4] that a slight decrease of ub, due to the widening band produces a slight increase 
in the probe velocity (blue curve). 

During experimentation, the granular medium is clearly seen to gain considerable freedom 
near the top surface and the beads closest to the faces of the probe tend to move away from 
the probe, reaching zones of the granular surface which hinder less the probe rotation. In 
doing so the beads leave the probe radially and, because of the asymmetry, contribute some 
net momentum to the probe and trigger rotation. The gained momentum is smaller for the 
faces parallel to the radii of the probe so the direction of rotation is that in which these 
faces move ahead. This is consistent with the observation that a completely immersed probe, 
no longer in contact with the top surface, ceases to rotate Fig. [2] (inset). In this respect 
the results are different from those reported in [l^ for the simulation of a Brownian motor 
driven by bacteria. 

Given the results obtained above, the following equation of motion for the probe may be 
written: 

lu = F - R - uboo + N (4) 
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where / is the moment of inertia of the rotating elements (the probe, its support and any 
dragged GM [l^), F is the net torque exerted by the GM on the probe, R is the frictional 
torque of the support, considered constant, ub^ represents the viscous resistance torque due 
to the granular p, [13], and is a noise term indicating fluctuations of F and is assumed 
to have zero time average, {N{t)) = 0. 

The main achievement of the experiment reported here has been the experimental obser- 
vation of the Brownian motor effect at a macroscopic, cm, scale. The observations are in 
general agreement with other works related to hydrodynamic viscosity [8|, ll7| and computer 
simulations [l^. A characteristic mechanical resonance is found in the range 30 — 40 Hz, 
similar to that found by numerical experiments [l^ which, when excited, triggers sponta- 
neous rotation of the immersed asymmetric probe. Excitation at higher frequency reduces 
the viscosity of the GM, but does not in itself trigger rotation. 
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